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We have presented theoretical and numerical studies on interactions of dispersive waves with 
second order solitons. Disperive wave with considerable intensity resonantly colliding with the 
second order solitons can lead to acceleration/decceleration of the later with subsequent central 
frequency shifts, but still well preserve the oscillating structure of the 2-solitons. The 2-soliton 
creates an effective periodical refractive index profile and thus reflected and transmitted dispersive 
waves generate new spectral bands with sharp peak structures. When resonant dispersive waves 
interact with two second order solitons bounding them, multiple scattering with multiple frequency 
conversion of the radiation occurs. Thus 2-solitons and radiation trapped in between them can 
produce effective solitonic cavity with ’’flat” or ’’concave mirrors” depending on the intensity of the 
input. 


Optical solitons attract attention of researches for 
decades due to the profound interest of their fundamen¬ 
tal physics and technological applications [T]. In addi¬ 
tion to the fundamental solitons, integrable models and 
physical media described by nearly-integrable equations 
[2] support high order A^-solitons, with N > 2, which 
are oscillating pulses periodically restoring their shape 
at distances that are multiples of the fundamental soliton 
period [3] . Experimental realization of 2- and 3-order op¬ 
tical solitons dates back to 1983 [4]. Initial narrowing of 
higher-order solitons was observed for soliton’s orders up 
to 13 [5]. Higher-order solitons were also demonstrated 
in the cavity of a mode-locked dye laser [6]. Strongly 
oscillating higher-order solitons find natural applications 
for the pulse compression HI ig and frequency conver¬ 
sion [9] , while the breakup of high order solitons initiates 
the extremely important process of the super continuum 
generation m- 

Propagation of optical pulse close to zero dispersion 
wavelength makes important consideration of high order 
dispersion terms and in particular the third order dis¬ 
persion (TOD) term. Second order solitons (2-solitons) 
were found to be practically robust under the influence of 
moderate TOD m , preserving their periodically oscillat¬ 
ing nature. Very recently resonant radiation of 2-solitons 
was investigated M revealing fascinating structure of 
radiation band consisting of multiple frequency peaks. 
While emitting the mentioned above radiation 2-soliton 
remained its fundamental features for many periods. In¬ 
terestingly, oscillating dissipative solitons featuring sim¬ 
ilar radiation behavior were reported too [lansi. 

In is relevant to mention that extensive studies were 
performed considering scattering of weak DWs on fun¬ 
damental solitons, with solitons properties remaining in¬ 
variant [HHSO]- On the other hand, it was also demon¬ 
strated that strong resonant collision of DWs with soli¬ 


tons can accelerate or decelerate solitons. These acceler¬ 
ations manifest themselves in frequency domain by up¬ 
shifts or downshifts of central frequencies of the solitons 
[2TII24] . Generation of new bands of frequencies due to 
the scattering of DWs from solitons can be viewed as an 
alternative technique for generation of broad and coher¬ 
ent supercontinuum [23] without the high order soliton 
fission. 

Taking one step further a very interesting scenario with 
DWs interacting with two co-propagating solitons 
ng. DWs temporally located in between the solitons in¬ 
teract with one soliton after the other bouncing of them. 
This process can repeat itself many times resulting in 
creation of solitonic cavity behaving like two mirrors for 
the DWs. In case of DWs with considerable intensity, 
the two solitons get attracted one towards the other, re¬ 
sulting in collision or even fusion [28] . This phenomenon 
was proven to be responsible for the appearance of mul¬ 
tiple soliton knot patterns [2TJ [26[ |29] during the com¬ 
plex super continuum generation process. Such a solitonic 
cavity was very recently successfully realized in experi¬ 
mental conditions Also a ’’convex mirror” cavity 

with mutually repulsing via DW two dark solitons was 
demonstrated [30] . 

In the present work we aim to study interactions of 
2-solitons with external dispersive waves (DWs), demon¬ 
strating a possibility to manipulate propagation trajecto¬ 
ries of the 2-solitons, to create broad continuum with rich 
spectral properties and also to create solitonic cavities 
consisting of two 2-solitons and DWs bouncing between 
them. 

Dynamics of optical pulse in the vicinity of zero dis¬ 
persion wavelength (ZDW) is governed by nonlinear 
Schrodinger equation with the inclusion of TOD. 

idzU - 1/32^4 M 71'wp'w = 0 (1) 
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where (32 and /^s designate second and third order dis¬ 
persion terms, while 7 represents the nonlinear parame¬ 
ter. Considering standard fiber with Xzdw = 1311 nm 
and launching a soliton at central wavelength of Xsoi = 
1470 nm, we will work with the following values of dis¬ 
persion parameters: ^2 = —14.2 ps^/km and Ps = 
0.087 ps^/km. The nonlinear coefficient is taken 7 = 
2 W“^km“^ 

The injected light consist of two pulses launched with 
different frequencies namely soliton and DW 

Uo{t) = Usol{t) + UDw{t) 

The 2-soliton has the initial form 

Usoi{t) = 2v^sech(t/To) 

with temporal width of Tq = 62.5 fs and the peak power 
of corresponding fundamental soliton Pq = 1817 W. DW 
is given by 

UDw{t) = A]jw sech((t — ti)/Ti) ex.p{—iSuj{t — ti)) 

Here Adw represent the amplitude of DW, ti the tem¬ 
poral delay of the DW with the respect to the soliton, Ti 
is the temporal width of dispersive wave, and Suj is the 
angular frequency deviation of the DW from the soliton 
given by Su = 27rc(A7 - Koi)- 

The result of collision of a relatively wide DW (Ti = 
2 OT 0 ) centered at 1130 nm and having a very small am¬ 
plitude of A‘jyy^ = 4.54 W with the 2 -sohton is presented 
in Fig. [^a, b). Fig. [^a) demonstrates dynamics in 
temporal domain and it reveals that the DW is partially 
scattered and partially transmitted through the 2 -soliton 
that creates an effective periodic refractive index poten¬ 
tial. In general one can archive much higher percentage 
of scattering by optimizing the incident wavelength of 
the DW within the resonance band, but here we aim to 
inspect also the transmitted part of the DW. Looking at 
the spectral domain representation of the dynamics in 
Fig. [^b) we can see several new frequency bands aris¬ 
ing in a course of evolution. The leftmost band below 
1100 nm consisting of multiple peaks is the resonant or 
Cherenkov radiation m emitted form the 2 -soliton. Its 
central spectral position as well as the separation between 
the frequency peaks is described by [14]: 

- Vg'^Suj = bp + 

2d ^2 Zq 

Here duj = 2ttc{X~^ — is frequency detuning of the 
resonant wave from the soliton frequency, Vg is the soliton 
velocity, Zq is a period of the oscillations of 2 -soliton and 
N is an integer number. 

The region of transmitted and reflected radiation at 
the output of the fiber is shown in Fig. [^c). The res¬ 
onance frequencies for the radiation generated by FWM 
mixing between incident DW and 2 -sohton were derived 


by using the perturbation theory similar to the analysis 
for the fundamental soliton in m- In contrast with the 
fundamental soliton case, the resonant radiation becomes 
polychromatic, with frequency detunings for incident and 
scattered fields duuinc, sc = - A“^^J being tied 

by the following relation 

( 2 ) 

We solved the <§ graphically in order to obtain the spec¬ 
tral positions of the transmitted and reflected peaks and 
presented it in Fig. 0d). For example the intersection 
of the blue curve with the horizontal magenta line rep¬ 
resents the scattered and transmitted DWs of the zero 
order (A^ = 0 ), that corresponds to the strongest re¬ 
flected peak at about 1223 nm and strongest transmitted 
peak at 1130 nm. The intersection of the blue curve with 
the horizontal green line stand for N = —1 and so on. 
Comparison between these results and those obtained by 
direct numerical simulations (Fig. [^c)) is visualized by 
vertical dashed lines. We can find an excellent agreement 
between the two approaches with a discrepancy of only 
few nanometers. 

Using more narrow incident DWs with the temporal 
width equal to that of the soliton (TI = TO) we have 
produced systematic simulations to provide the depen¬ 
dence of the scattered radiation central wavelength vs. 
the incident central wavelength as shown in Fig. Solid 
curve is the plot of the analytical relation Xsc(Xinc) that 
is obtained from graphical solution of ([^ for the funda¬ 
mental harmonic N = 0. Shifting ti and thus modifying 
the place of collision with the respect to soliton’s shape at 
the moment of the collision did not significantly influence 
the results. 

Now we will inspect what happens when DW with 
more considerable intensity collides with the 2-soliton. In 
this case the velocity of the soliton will be modified and 
we can no longer rely on analytical approximations, but 
to study only results of numerical simulations. We collide 
DW with the 2 -sohton as in the previous chapter, but this 
time we increase its input intensity to = 454.4 W, 

that is still several times smaller that the input peak in¬ 
tensity of the 2 -sohton participating in the collision and 
we will take the width of the DW equal to that of the 
soliton (Ti = To). A typical result of collision of strong 
DW with the 2 -soliton still resisting to the splitting is 
presented in the Fig. 3. The soliton’s trajectory is con¬ 
siderably shifted, consequently it central wavelength is 
downshifted from 1470 nm to 1463 nm. We can also 
observe in Fig. 3 (lower panel) generation of significant 
new light spectral bands in region between the soliton’s 
central wavelength and DW. Increasing the intensity of 
the incident DW much higher we will be observing the 
splitting of the 2 -sohton. 

Even more interesting interaction scenario occurs when 
DW bounces off two co-propagating 2 -sohtons. We con- 
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FIG. 3: (Color online) Interaction of a strong DW with 2- 
soliton. Dynamics in (a) temporal domain showing (b) 

Spectral domain dynamics representation in logarithmic scale. 


sider an input consisting of two well separated 2-solitons 
and DW launched in between them such as 


FIG. 1: (Color online) Scattering of a wide DW on 2-soliton. 
Dynamics in (a) temporal domain with \u\^'^ shown instead 
of the conventional intensity in order to observe better less 
intense regions, (b) Spectral domain dynamics presented in 
logarithmic scale, (c) The spectral density at the end out¬ 
put of the fiber with the region of reflected and transmitted 
DW zoomed (from 1100 nm to 1255 nm). (d) is the graph¬ 
ical solution to § predicting the position of the resonance 
frequencies. Blue curve is the left hand side of the equa¬ 
tion; horizontal lines correspond to the right hand’s side for 
N = — 1, 0,.., 3 starting from the bottom. 



FIG. 2: (Color online) Central wavelength of the scattered 
radiation as a function of incident. Solid line is the curve 
predicted by § for fundamental {N = 0) harmonic. Square 
markers are found by numerical modeling. 


uo{t) = 2vF)Sech((t + ti)/ro) 

+ A]jw sech(t/To) ex.p{—i6ujt) 

+ 2v^sech((t - ti)/To) 

If the wavelength of DW is located well inside the res¬ 
onance band we will be able to observe its multiple re¬ 
flections from the 2-solitons If the intensity of the DW 
is not very high the trajectories and the spectral char¬ 
acteristics of the two solitons are almost not modified 
and reflections occur pretty periodically to certain wave¬ 
length regions. For example the DW initially launched at 
1140 nm with = 18.17 W is being reflected several 
times to the region around 1200 nm and then back to 
its original spectral location (Fig. [^a, b)). Interestingly 
the spectrum of the two solitons manifests multiple pe¬ 
riodically alternating light and dark regions spanned for 
this particular example from 1250 nm to 1750 nm. Such 
an illumination spectrum that can be further controlled 
by the parameters of the participating solitons can be 
suggested for photonics applications requiring particular 
spectral design. 

As mentioned above 2-solitons feature high degree of 
robustness when colliding with DWs and thus we can 
build a closed cavity consisting of two 2-solitons acting 
like concave mirrors without observing their splitting. An 
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FIG. 4: (Color online) Solitonic cavity created by two 2- 
solitons and DW bouncing between them, (a) Temporal do¬ 
main representation with shown, (b) Spectral domain 

representation with output spectrum profile shown above the 
panel. 


example of such cavity is demonstrated in Fig. ia, b). 
Here we used the same type of input light but increased 
an intensity of the incident DW to = 90.85 W 

leaving the rest of the parameters unchanged. In this 
example the central wavelength of the left and right 2- 
solitons get upshifted and downshifted respectively. Also 
the periodicity of illuminated regions such as in Fig. 
gets strongly distorted (Fig. [^b)). 
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FIG. 5: (Color online) Solitonic ’’concave mirror” cavity cre¬ 
ated by two 2-solitons and stronger DW bouncing between 
them, (a) Temporal domain representation with shown, 
(b) Spectral domain representation. 

To conclude, scattering of DWs on high order solitons 
was presented. The structure of the 2-soliton provides 


an effective periodical potential leading to polychromatic 
scattering and transmission for the incident DWs. In 
spectral domain it manifest itself by the generation of 
frequency comb like novel bands. Position of the peaks of 
the transmitted and reflected bands is predicted by excel¬ 
lently overlapping analytical and numerical approaches. 
Interaction of DWs with two 2-solitons bounding them 
and creation of solitonic cavities was also demonstrated. 
The presented nonlinear dynamics besides the fundamen¬ 
tal interest suggest prospective photonic devices schemes 
for frequency conversion and broadband light generation 
with pre-designed spectrum. 
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